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The prevailing view at present is that postsynaptic expression of the classical 
NMDA receptor-dependent long-term potentiation relies on an increase in 
the numbers of local AMPA receptors (AMPARs). This is thought to parallel 
an expansion of postsynaptic cell specializations, for instance dendritic spine 
heads, which accommodate synaptic receptor proteins. However, glutamate 
released into the synaptic cleft can normally activate only a hotspot of 
low-affinity AMPARs that occur in the vicinity of the release site. How the enlar- 
gement of the AMPAR pool is causally related to the potentiated AMPAR 
current remains therefore poorly understood. To understand possible scenarios 
of postsynaptic potentiation, here we explore a detailed Monte Carlo model of 
the typical small excitatory synapse. Simulations suggest that approximately 
50% increase in the synaptic AMPAR current could be provided by expanding 
the existing AMPAR pool at the expense of 100-200% new AMPARs added at 
the same packing density. Alternatively, reducing the inter-receptor distances 
by only 30-35% could achieve a similar level of current potentiation without 
any changes in the receptor numbers. The NMDA receptor current also appears 
sensitive to the NMDA receptor crowding. Our observations provide a quanti- 
tative framework for understanding the 'resource-efficient' ways to enact 
use-dependent changes in the architecture of central synapses. 



Keywords: 

AMPA receptor, NMDA receptor, long-term 
potentiation 

Author for correspondence: 

Dmitri A. Rusakov 
e-mail: d.rusakov@ucl.ac.uk 



1. Introduction 

Cellular mechanisms of use-dependent synaptic plasticity remain a subject of 
intense investigation, simply because they hold a promise to unveil the neuro- 
biological basis of learning and memory formation in the brain. The molecular 
machinery underlying the classical plasticity paradigm, long-term potentiation 
(LTP) of AMPA receptor (AMPAR)-mediated synaptic transmission, has been at 
the centre stage of neuroscience research for decades [1-3]. The currently 
prevailing view is that the expression of classical NMDA receptor (NMDAR)- 
dependent LTP relies on the increased postsynaptic AMPAR current, although 
a boost in presynaptic release probability has been found important at least 
in some physiological scenarios [4-6]. The mechanism providing the LTP- 
associated AMPAR current increase is thought to involve activity-induced 
insertion of additional synaptic AMPARs through either endocytosis or lateral 
membrane trafficking, or both [7-12]. This paradigm could, in principle, 
explain a variety of physiological phenomena associated with a use-dependent 
functional and structural plasticity of excitatory synaptic connections [13]. 

At the same time, various theoretical models of small excitatory synapses 
seem to converge on a prediction that glutamate released into the synaptic 
cleft from a synaptic vesicle activates in most cases only a relatively small hot- 
spot of local AMPARs [14-19]. Consequently, it has long been understood that 
the pattern of AMPARs within the synaptic cleft could significantly affect 
the synaptic current amplitude. However, whether morphologically plausible 
changes in the receptor distribution inside the cleft could actually produce 
synaptic potentiation consistent with the experimental LTP has remained 
unclear. Similarly, the question arises whether the insertion of additional 
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Figure 1. Simulations of glutamate release, diffusion and receptor activation in the environment of small central synapses, (a) Model schematic; see §2tf for notation 
and [20,22] for algorithm details and experimental validation, (b) An example of AMPAR current simulations in response to release of 2700 glutamate molecules 
(postsynaptic membrane potential V m = -80 mV). Grey staggered traces, 10 example simulation runs depicting stochastic activation of individual AMPARs by 
individual glutamate molecules, as further detailed in [18]; black trace, an average trace of 40 simulation runs, with stochastic receptor activation; red trace, simu- 
lation outcome with AMPARs activated by the average glutamate concentration calculated from the number of glutamate molecules in the vicinity (to reduce 
computation time). 



AMPARs simply by expanding the postsynaptic density 
(PSD), as some LTP scenarios appear to suggest, is a naturally 
efficient way to increase synaptic strength. 

Obtaining reliable experimental evidence for an LTP- 
associated increase in the local AMPAR numbers or their density 
on the nanoscale has not yet been technically feasible. Where, 
within the PSD, and how many additional AMPARs are 
normally required to explain robust synaptic potentiation 
remains therefore poorly understood. To address these questions 
in a quantitative manner, here we explore the relationship 
between the arrangement of AMPARs (and NMDARs) and 
synaptic efficacy at common excitatory synapses (exemplified 
by the CA3-CA1 connection in the hippocampus) using a 
computational Monte Carlo model. We adapt the modelling 
approach which has been tested and validated extensively 
against the experimental data including sub-millisecond prob- 
ing of native AMPARs and NMDARs in outside-out, somatic 
and dendritic, membrane patches [18,20-23]. In our simulations, 
we systematically explore changes in the density and numbers of 
synaptic AMPARs (and NMDARs) and examine the ensuing 
effect on the amplitude of synaptic currents. We have to stress 
here that this study focuses on a widely held, archetypal 
LTP paradigm which is confined to an individual synapse and 
has been routinely explored in various experimental circum- 
stances, from electron microscopy to dendritic spine imaging 
studies. The LTP mechanisms which involve changes in the 
composition of synaptic population, for instance conversion of 
silent synapses, are outside the scope of the present analysis. 
In assessing the reliability of our theoretical predictions, we rou- 
tinely test whether such predictions remain robust over a 
physiological range of synaptic parameters that are difficult to 
access empirically. 



2. Material and methods 

(a) Monte Carlo model: synaptic architecture and 
environment 

We adapted the modelling methods which have been tested and 
validated previously with respect to the architecture and physi- 
ology of CA3-CA1 synapses [18,20,22,24], Computations were 



carried out using an in-house 64-node PC cluster optimized for 
parallel computing [18] (algorithms provided by Sitrus LLC, 
Boston, MA, USA). In the model, the synaptic apposition between 
adjacent pre- and postsynaptic membranes was represented by a 
flat cylindrical cleft with radius R ranging from 180 to 240 nm 
(figure la). This included the PSD (the membrane area which 
contains randomly scattered synaptic AMPARs and/or 
NMDARs and is opposed by the presynaptic active zone 
where release occurs) with radius r a ranging from 60 to 200 nm 
(figure la). The diffusion coefficient of glutamate inside the 
cleft was set at 0.33 |xm 2 ms -1 , and 0.4 |mm 2 ms -1 outside the 
cleft, in accord with previously published estimates [18,22,25]. 
The synaptic cleft height 8 = 19 nm was set to reflect the average 
distance between the electron density maxima of the pre- and 
postsynaptic membranes in electron micrographs of CA1 
synapses [20,26]. Within the postsynaptic active zone, the total 
number N of AMPARs (channel conductance 12.5 pS) or 
NMDARs (25 pS) varied between 15 and 110 or 2 and 30, 
respectively, consistent with quantitative immuno-electron 
microscopy [27,28] and with the range estimated by optical quan- 
tal analyses [29] at these synapses. Brownian movements of 
individual glutamate molecules and receptor state transitions fol- 
lowing activation were computed with a time step of 0.1 |uls; 
further reduction of the time step by an order of magnitude 
improved computation accuracy by less than 1%; all kinetic 
receptor states were recorded and stored. 

(b) Monte Carlo model: release of glutamate 

A total of 2700 glutamate molecules were released into the cleft 
quasi-instantaneously, according to the value estimated in a 
recent experimental study of CA3-CA1 synapses [22]. In most 
cases, glutamate was released at the cleft centre, in accord with 
a traditional modelling approach that routinely accounts for the 
problem's rotational symmetry. However, we also explored the 
scenario in which the release site was distributed randomly 
within the active zone. In these latter simulations, the average 
amplitude of the AMPAR or NMDAR currents (and its standard 
deviation) was calculated for 20 consecutive simulated release 
events, as indicated. 

(c) Monte Carlo model: receptor activation 

The working algorithm, which has been detailed earlier [22,24,30], 
calculates receptor kinetics from average concentrations of local 
glutamate represented by individually tracked molecules in 



the vicinity of the receptor. This hybrid algorithm requires 
much less computational resources compared with the Monte- 
Carlo-throughout method which computes explicit stochastic 
behaviour of individual receptors activated by individual gluta- 
mate molecules [18] while giving a virtually identical average 
outcome (figure lb). In brief, at each time step (dt = 0.1 |xs), the 
model updated the coordinates of all individual glutamate mol- 
ecules that follow Brownian movement. Next, it calculated the 
concentration profile of glutamate C(r,t) in the cleft. In conditions 
of approximate rotational symmetry, this corresponded to 
C(r,t) = N 8 (27rr8Ar)~ 1 , where N 8 stands for the number of gluta- 
mate molecules occurring at time point t inside the flat 
cylindrical ring of height 8, width Ar and radius r. The average 
occurrence (concentration) of open receptors [0](r) within the 
active zone (r<r a ) was then calculated at the same time point 
from the multi-stage AMPAR kinetics [31] or NMDAR kinetics 
[32] adjusted for 33-35°C using the immediate history of the recep- 
tor states and C(r,t). These calculations gave the total synaptic 
current (V m =-80mV for AMPARs, no Mg 2+ block for 
NMDARs) integrated over the PSD area. The duty cycle was 
repeated systematically throughout the model run. 

3. Results 

(a) Potentiation requires a disproportionately large 
number of AMPARs added to the PSD at the 
same density 

We first asked how rapidly the synaptic AMPAR current 
increases with new AMPARs added to the PSD (with the 
PSD expanding in area as new AMPARs are added) assum- 
ing that the average surface density ('packing') of AMPARs 
remains unchanged. Simulation results indicated that, for a 
variety of synaptic cleft sizes (R), the relationship between 
the number of added AMPARs and the synaptic current is 
relatively weak: to increase the current by approximately 
50%, the receptor numbers have to be boosted by 
100-200% (figure 2a). This relationship was qualitatively 
similar for various cleft dimensions and across the range of 
AMPAR pool (PSD) sizes. The simple reason for such a dis- 
proportionally weak effect is that AMPARs added to the 
PSD periphery are much less sensitive to glutamate released 
near the PSD centre compared with AMPARs expressed 
closer to the release site. These simulations also indicated 
that, in principle, expanding the synaptic cleft (R) on its 
own could also boost the synaptic current, with the effect 
growing stronger at higher AMPAR densities. Again, this 
relationship is due to the simple fact that the two-dimensional 
cleft corresponds to greater (i.e. stronger and longer) 
local transients of released glutamate compared to the 
three-dimensional extracellular space surrounding it. 

(b) Modest crowding of synaptic AMPARs on the 
nanoscale can boost synaptic current without any 
new receptors 

We next changed the scenario and asked whether the contraction 
of the synaptic AMPAR pool, i.e. a reduction in the average 
distance between receptors without any changes in the receptor 
numbers, would produce detectable changes in the total 
AMPAR current. Perhaps surprisingly, our data predicted that 
approximately 50% boost in the total current requires a relatively 
modest decrease in the inter-receptor distance (30-35%; 



figure 2b). Again, this effect remained robust for different 
synaptic cleft dimensions and over a wide range of AMPAR 
pool (PSD) sizes that reflect the average nearest-neighbour dis- 
tances between individual AMPARs. Alternatively, inserting 
new receptors to the PSD without changing the PSD size, i.e. 
increasing both AMPAR number and their surface density, had 
a combined boosting effect on the AMPAR current (figure 2c). 

(c) Varying the intra-cleft location of glutamate release 

Although in the above simulations the glutamate release site 
was placed at the cleft centre, this may not necessarily be the 
case in reality: CA3-CA1 synapses display a range of locations 
for docked (ready-releasable) synaptic vesicles within the pre- 
synaptic active zone [34]. Varying the release site location in a 
synapse model was previously demonstrated to reduce the 
average synaptic current while substantially increasing its 
amplitude variability [16]. To assess the effect of this variation 
on our conclusions, we explored the characteristic scenario in 
which the typical synapse (r a = 110 nm, R = 200 nm, N = 50 
AMPARs) underwent approximately 50% potentiation either 
by adding AMPARs at the same density or by increasing 
AMPAR density without changing their number. We therefore 
first looked at this scenario with the centre-fixed glutamate 
release site and then repeated simulations for exactly the 
same receptor arrangements and the same synaptic architec- 
ture but with the release site being randomly distributed, 
trial-to-trial, over the presynaptic active zone which opposes 
the PSD. The comparison showed that the stochastic occur- 
rence of the release site generally decreased the AMPAR 
current amplitude while substantially increasing its variability 
(figure 2d), which was consistent with previous estimates [16]. 
At the same time, however, the effect of AMPAR pool expan- 
sion or clustering on the synaptic current potentiation (approx. 
50% LTP) was generally indistinguishable from that observed 
in the centre-fixed release scenario (figure 2d). 

(d) A certain combination of modest changes in 
receptor density and numbers alters synaptic 
efficacy 'most efficiently' 

To evaluate and compare relative contributions of various 
changes in receptor arrangement (and synaptic configur- 
ation) leading to potentiation or depression of the AMPAR 
current, we explored the synaptic parameter space more 
systematically. The resulting parametric map (figure 2e) 
helps us to understand the pattern of synaptic current altera- 
tions pertinent to changes in the AMPAR pool (PSD) size or 
in the total AMPAR number (and thus the average inter- 
receptor distance), or in both parameters simultaneously. 
For instance, these data illustrate quantitatively to what 
degree the synaptic AMPAR current is sensitive to a varia- 
tion in the AMPAR numbers at a constant surface density 
(figure 2e, curved dotted arrows marked Density = const), or 
to receptor clustering with no extra receptors (figure 2e, hori- 
zontal dotted arrows marked N = const). They also indicate 
that a certain type of synaptic receptor rearrangement can 
alter the synaptic response in a most 'economical' way, i.e. 
by implementing the smallest, in relative terms, structural 
alteration of the synapse. This type of change corresponds 
to the steepest gradient for each data point which reflects a 
given synaptic configuration (figure 2e; exemplified by the 
arrow marked fastest). 
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Figure 2. Potentiation of the synaptic AMPAR current can be achieved either by adding many more AMPARs to the pool, or by relatively modest clustering (crowding) of 
existing AMPARs, or both, (a) (i) Diagrams (not to scale) and traces: example simulated AMPAR currents before (black dots, trace) and after (grey) adding approximately 
130% extra receptors to existing 50 AMPARs by expanding the receptor pool (PSD) from r a = 110 nm to r a = 160 nm (unchanged receptor density, synaptic cleft radius 
R = 200 nm). (ii) The relationship between the number of local synaptic AMPARs (or the PSD size) and the total AMPAR current under the unchanged surface density of 
AMPARs, for four synaptic cleft sizes /?, as indicated; dots, results of individual simulation experiments; coloured dashed lines, linear regression; grey vertical dotted line 
(here and thereafter) indicates the experimental average size of the PSD at hippocampal CA3-CA1 synapses [26,33]; two dotted arrows illustrate an arbitrarily selected 
example of changes that corresponds to an approximately 50% potentiation of the AMPAR current, (b) (i) Diagrams (not to scale) and traces: example simulated AMPAR 
currents before (black trace) and after (grey) shrinking the pool of N = 50 AMPARs from r a = 110 nm to r a = 70 nm (PSD radius; R = 200 nm), with no changes in N. 
The relationship between the average inter-receptor nearest-neighbour distance (or the AMPAR pool size) and the total AMPAR current under the unchanged total number 
of AMPARs; coloured dashed lines, quadratic regression; other notations as in (a), (c) An example of the relationship between the synaptic cleft size (/?), the total number of 
synaptic AMPARs and the total AMPAR current with the PSD (AMPAR pool) size remaining unchanged at r a = 110 nm. (d) Testing the impact of glutamate release site 
variability. Red circles (mean + s.d., n = 20 trials), centre-fixed release site data (as in a-c) for the typical synapse in control conditions (cntrl, r a = 110 nm, 
R = 200 nm, N = 50 AMPARs) and in conditions that correspond to approximately 50% AMPAR current potentiation, either through insertion of additional AMPAR 
to the periphery (LTP insert, r a = 180 nm, R = 200 nm, N = 134 AMPARs) or by shrinking the PSD without any other changes (LTP clust, r a = 75 nm). Green circles 
(mean + s.d., n = 20 trials), data obtained under evenly random variation of the release site location over the active zone (opposite to the PSD); other synaptic and 
receptor parameters are the same as in control conditions (cntrl). (e) Parametric map depicting the relationship between the AMPAR pool size, receptor density and total 
AMPAR current (false colour scale). Circle, a reference configuration corresponding to the experimental average (r a = 110 nm, R = 200 nm, N = 50 AMPARs); horizontal 
dotted arrows show how a change in the receptor density alone affects AMPAR current; curved dotted arrow shows how changes in the receptor numbers (unchanged 
density) affect the AMPAR current; red dotted arrows depict the direction of the fastest change. 



(e) Synaptic responses are also sensitive to NMDAR 
clustering 

The estimates above suggest that simply expanding the PSD by 
adding new AMPARs is a considerably less 'efficient' way to 
boost the synaptic strength when compared with increases in 



the AMPAR density. Because many types of excitatory synapses 
express another common ionotropic glutamate receptor, 
NMDAR, it was important to understand what changes in 
the NMDAR expression and pattern could potentiate the synap- 
tic NMDAR current. We therefore carried out a theoretical 
exploration of NMDAR activation by glutamate in conditions 
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Figure 3. Potentiation of the synaptic NMDAR current can be achieved either by adding more NMDARs to the pool, or by clustering (crowding) of existing NMDARs, 
or both, (a) (i) Traces: example simulated NMDAR currents before (black) and after (grey) adding 12 new receptors (120% extra receptors) to existing 10 NMDARs by 
expanding the receptor pool from r a = 110 nm to r a = 165 nm (NMDAR density unchanged, R = 200 nm). (ii) The relationship between the number of local 
synaptic NMDARs (or the NMDAR pool size) and the total NMDAR current, with the unchanged NMDAR density. Two dotted arrows illustrate an arbitrarily selected 
example of changes that corresponds to approximately 50% potentiation of the NMDAR current. Other notations are as in figure 2a. (b) (i) Traces: example simulated 
NMDAR currents before (black) and after (grey) shrinking the pool of N = 10 NMDARs from r a = 110 nm to r a = 73 nm (R = 200 nm), with no changes in N. The 
relationship between the average inter-receptor nearest-neighbour distance and the total NMDAR current under the unchanged total number of NMDARs; coloured 
dashed lines, quadratic regression. Other notations are as in (a), (c) An example of the relationship between the synaptic cleft size (/?), the total number of synaptic 
NMDARs and the total NMDAR current with the PSD (NMDAR pool) size remaining unchanged at r a = 110 nm. (d) Testing the impact of glutamate release 
site variability. Red circles (mean + s.d., n = 20 trials), centre-fixed release site data (as in a-c) for the typical synapse in control conditions (cntrl, r a = 
110 nm, R = 200 nm, N = 10 NMDARs) and in conditions that correspond to approximately 50% NMDAR current potentiation, either through insertion of additional 
NMDARs to the periphery (LTP insert, r a = 190 nm, R = 200 nm, N = 30 NMDARs) or by shrinking the PSD without any other changes (LTP dust, r a = 80 nm). Green 
circles (mean + s.d., n = 20 trials), data obtained under evenly random variation of the release site location over the active zone (opposite to the PSD); other synaptic 
parameters are the same as in control conditions (cntrl). (e) Parametric map depicting the relationship between the NMDAR pool (PSD) size, NMDAR density and total 
NMDAR current (colour scale); circle, a reference configuration corresponding to the experimental average (r a = 110 nm, R = 200 nm, N = 10 NMDARs); horizontal 
dotted arrows show how a change in the receptor density alone affects NMDAR current; curved dotted arrow shows how changes in the receptor numbers (unchanged 
density) affect the NMDAR current; red dotted arrows depict the direction of the fastest change. 



similar to those we considered for AMPARs (figure 3). Although 
activation of higher affinity NMDARs is supposed to be much 
less sensitive to the distance from the glutamate release site com- 
pared with lower affinity AMPARs [18,35,36], our simulation 



results indicate that this distinction between AMPARs and 
NMDARs is less prominent inside the cleft. Indeed, the 
NMDAR current was also sensitive to both adding NMDARs 
to the PSD and to NMDAR clustering, although the difference 



between the two effects was somewhat less striking compared 
with that in the case of AMPARs (figure 3a- c). Similar to the 
case of AMPARs, varying the release site location reduced the 
average amplitude of NMDAR current while increasing its varia- 
bility (figure 3d). At the same time, the effect pertinent to LTP 
expression remained qualitatively similar (figure 3d). Again, a 
systematic exploration of synaptic parameters provided intuitive 
clues as to what type of changes in the NMDAR expression or 
pattern could produce the fastest, and thus arguably the most 
'efficient', alteration in the NMDAR current (figure 3e). 



4. Discussion 

Here, we have employed a Monte Carlo model of small central 
synapses to examine what type of changes in local AMPAR 
(and NMDAR) expression might explain the receptor current 
increase following the induction of LTP, assuming no changes 
in the amount of glutamate released. The results indicate that 
the addition of new AMPARs simply by expanding the PSD 
results in a disproportionally small boost of the AMPAR current. 
On average, a given fractional increase in the AMPAR current 
requires a three to four times larger increase in the relative 
number of AMPARs added this way. By contrast, a relatively 
modest (approx. 30%) reduction in the average inter-AMPAR 
distances within the PSD, and no additional AMPARs, could 
provide a robust (approx. 50%) rise in the AMPAR current 
amplitude. Our tests suggested that this conclusion was largely 
independent of whether the glutamate release site was fixed at 
the cleft centre or varied randomly over the active zone. This out- 
come was held over a wide range of poorly controlled synaptic 
parameters, suggesting that the conclusion was relatively 
robust. Importantly, we also found that for a given synaptic 
architecture, there is a certain combination of changes in 
AMPAR (and NMDAR) density and numbers that could provide 
the 'most economical', in terms of the relative structural change, 
path to the required level of synaptic potentiation or depression. 

Notwithstanding the lack at present of direct experimental 
evidence supporting the hypothesis of use-dependent AMPAR 
nanoscale clustering, there are several factors, in addition to 
being by far the most parsimonious explanation, which make 
this hypothesis intellectually attractive. Firstly, this mechanism 
is readily consistent with the fact that LTP reflects an increased 
number of activated postsynaptic AMPARs. Because AMPARs 
are normally far from saturation at the synapses under study 
[22,37-39], the expression of LTP might, in principle, reflect 
either activation of newly inserted AMPARs or a higher acti- 
vation level of the existing AMPARs, or perhaps both. 
Secondly, because AMPAR clustering does not require an 
increase in the total receptor numbers expressed within the 
PSD, the involvement of complex cellular machineries engaged 
in use-dependent AMPAR trafficking might not necessarily be 
critical throughout all stages of LTP expression. In this context, 
it is noteworthy that manifestation of AMPAR trafficking is 
usually detected minutes after the plasticity-inducing stimulus, 



whereas electrophysiological indicators point to a virtually 
immediate synaptic efficacy increase following LTP induction. 
Indeed, it is reasonable to think that structural rearrangement 
on the nanoscale is likely to occur more rapidly than changes 
associated with molecular transport over macroscopic dendritic 
compartments. Thirdly, the fact that some of the molecular cas- 
cades responsible for AMPAR trafficking have been found as a 
prerequisite for LTP expression as such might also reflect an 
important role of the same molecular machinery for nanoscale 
AMPAR clustering. 

Our results suggest that activation of high-affinity NMDARs 
is also sensitive to the clustering inside the synaptic cleft. It is 
tempting to speculate that this mechanism could, in principle, 
contribute to activity-dependent changes in NMDAR-mediated 
postsynaptic Ca 2+ entry. The latter in turn could affect local 
conditions for Ca 2+ -dependent plasticity. NMDARs are thought 
to undergo much slower macro- and microscopic movements 
compared with AMPARs, also depending on the receptor iso- 
form [40,41]. Their 'reluctance' to move suggests that small 
nanoscale rearrangement of NMDARs could play a significant 
role in plasticity changes involving NMDAR signalling. 

The features of synaptic receptor crowding on the nanoscale 
are presently beyond the resolution of optical microscopy or 
other real-time recording methods, making it difficult to test 
and explore the underlying machinery in a direct fashion. None- 
theless, the typical PSD of central synapses appears to have 
enough space to accommodate variable numbers of AMPARs 
and NMDARs at a range of surface densities [42]. This suggests 
that use-dependent changes in receptor crowding within the 
PSD are not implausible. Interestingly, the key PSD constituent, 
the PSD-95 protein, appears to give rise to scaffolding-like 
filamentous links that space synaptic AMPARs and NMDARs 
at relatively regular 20-30 nm intervals [43], which appear 
somewhat greater than the inter-receptor distances explored in 
our simulations. However, freeze-fracture replica labelling 
studies reveal smaller intervals between individual AMPARs 
in the PSD [44]. More importantly, our simulations quote the 
average nearest-neighbour distance assuming a uniformly 
random receptor distribution (Poisson point process), whereas 
this characteristic distance is substantially greater for a regular 
expression pattern, for instance a two-dimensional square lattice, 
with the same surface density [45]. The aim of our study was to 
explore several hypotheses (pertinent to the microscopic mech- 
anisms of LTP expression) that would appear physiologically 
plausible and yet consistent with the existing experimental 
evidence. It would require, however, a new generation of experi- 
mental techniques to test what actually controls the receptor 
expression pattern inside the synaptic cleft and whether this 
pattern indeed changes in the course of use-dependent plasticity. 



Funding statement. This work was supported by the Wellcome Trust Prin- 
cipal Fellowship (UK), European Research Council Advanced grant, 
Medical research Council (UK), Biology and Biotechnology Research 
Council (UK), BM1001 Cost Action (European Commission). 



to 
b 
<5 



vO 



References 

1. Bliss T, Lomo T. 1973 Long-lasting potentiation of 
synaptic transmission in the dentate area of the 
anaesthetized rabbit following stimulation of the 
perforant path. J. Physiol. Ill, 331-356. 



2. Malenka RC, Nicoll RA. 1999 Long-term 
potentiation: a decade of progress? Science 
285, 1870-1874. (doi:10.1126/science.285. 
5435.1870) 



3. Bliss T, Collingridge G. 1993 A synaptic model of 
memory: long-term potentiation in the 
hippocampus. Mure 361, 31-39. (doi:10.1038/ 
361031a0) 



4. Malgaroli A, Tsien RW. 1992 Glutamate-induced 
long-term potentiation of the frequency of 
miniature synaptic currents in cultured hippocampal 
neurons. Nature 357, 134-139. (doi:10.1038/ 
357134a0) 

5. Kullmann DM, Erdemli G, Asztely F. 1996 LTP of 
AMPA and NMDA receptor-mediated signals: 
evidence for presynaptic expression and 
extrasynaptic glutamate spill-over. Neuron 17, 
461 -474. (doi:10.1016/S0896-6273(00)80178-6) 

6. Emptage NJ, Reid CA, Fine A, Bliss TV. 2003 Optical 
quantal analysis reveals a presynaptic component of 
LTP at hippocampal Schaffer-associational synapses. 
Neuron 38, 797-804. (doi:10.1016/S0896- 
6273(03)00325-8) 

7. Shi SH, Hayashi Y, Petralia RS, Zaman SH, Wenthold 
RJ, Svoboda K, Malinow R. 1999 Rapid spine 
delivery and redistribution of AMPA receptors after 
synaptic NMDA receptor activation. Science 284, 

1 81 1 - 1 81 6. (doi: 1 0. 1 1 26/science.284.5421 .1811) 

8. Noel J, Ralph GS, Pickard L, Williams J, Molnar E, 
Uney JB, Collingridge GL, Henley JM. 1999 Surface 
expression of AMPA receptors in hippocampal 
neurons is regulated by an NSF-dependent 
mechanism. Neuron 23, 365-376. (doi:10.1016/ 
S0896-6273 (00)80786-2) 

9. Borgdorff AJ, Choquet D. 2002 Regulation of AMPA 
receptor lateral movements. Nature 417, 649-653. 
(doi : 1 0. 1 038/nat u re00780) 

10. Park M, Penick EC, Edwards JG, Kauer JA, Ehlers MD. 
2004 Recycling endosomes supply AMPA receptors 
for LTP. Science 305, 1972-1975. (doi:10.1 126/ 
science.1 102026) 

11. Plant K, Pelkey KA, Bortolotto ZA, Morita D, 
Terashima A, McBain CJ, Collingridge GL, Isaac JT. 
2006 Transient incorporation of native GluR2-lacking 
AMPA receptors during hippocampal long-term 
potentiation. Nat. Neurosci. 9, 602-604. (doi:10. 
1038/nn1678) 

12. Makino H, Malinow R. 2009 AMPA receptor 
incorporation into synapses during LTP: the role 
of lateral movement and exocytosis. Neuron 64, 
381-390. (doi:10.1016/j.neuron.2009.08.035) 

13. Czondor K, Mondin M, Garcia M, Heine M, 
Frischknecht R, Choquet D, Sibarita JB, Thoumine 
OR. 2012 Unified quantitative model of AMPA 
receptor trafficking at synapses. Proc. Natl Acad. Sci. 
USA 109, 3522-3527. (doi:10.1073/pnas. 
1109818109) 

14. Uteshev VV, Pennefather PS. 1996 A mathematical 
description of miniature postsynaptic current 
generation at central nervous system synapses. 
Biophys. J. 71, 1256-1266. (doi:10.1016/S0006- 
3495(96)79325-6) 

15. Rusakov DA. 2001 The role of perisynaptic glial 
sheaths in glutamate spillover and extracellular 
Ca 2+ depletion. Biophys. J. 81, 1947-1959. 
(doi : 1 0. 1 01 6/S0006-3495 (01 )75846-8) 

16. Franks KM, Stevens CF, Sejnowski TJ. 2003 
Independent sources of quantal variability at 
single glutamatergic synapses. J. Neurosci. 23, 
3186-3195. 



17. Lisman JE, Raghavachari S, Tsien RW. 2007 The 
sequence of events that underlie quantal 
transmission at central glutamatergic synapses. Nat. 
Rev. Neurosci. 8, 597-609. (doi:10.1038/nrn2191) 

18. Zheng K, Scimemi A, Rusakov DA. 2008 Receptor 
actions of synaptically released glutamate: the role 
of transporters on the scale from nanometers to 
microns. Biophys. J. 95, 4584-4596. (doi:10.1529/ 
biophysj.108.129874) 

19. Rusakov DA, Savtchenko LP, Zheng K, Henley JM. 
2011 Shaping the synaptic signal: molecular 
mobility inside and outside the cleft. Trends 
Neurosci. 34, 359-369. (doi:10.1016/j.tins. 
2011.03.002) 

20. Savtchenko LP, Rusakov DA. 2007 The optimal 
height of the synaptic cleft. Proc. Natl Acad. Sci. USA 
104, 1823-1828. (doi:10.1073/pnas.0606636104) 

21. Sylantyev S, Savtchenko LP, Ermolyuk Y, Michaluk P, 
Rusakov DA. 2013 Spike-driven glutamate 
electrodiffusion triggers synaptic potentiation via a 
Homer-dependent mGluR-NMDAR link. Neuron 77, 
528-541. (doi:10.1016/j.neuron.2012.11.026) 

22. Savtchenko LP, Sylantyev S, Rusakov DA. 2013 
Central synapses release a resource-efficient amount 
of glutamate. Nat. Neurosci. 16, 10-12. (doi:10. 
1038/nn.3285) 

23. Sylantyev S, Rusakov DA. 2013 Sub-millisecond 
ligand probing of cell receptors with multiple 
solution exchange. Nat. Protoc. 8, 1299-1306. 
(doi:10.1038/nprot.2013.075) 

24. Sylantyev S, Savtchenko LP, Niu YP, Ivanov Al, 
Jensen TP, Kullmann DM, Xiao MY, Rusakov DA. 
2008 Electric fields due to synaptic currents sharpen 
excitatory transmission. Science 319, 1845-1849. 
(doi:10.1126/science.1 154330) 

25. Nielsen TA, DiGregorio DA, Silver RA. 2004 
Modulation of glutamate mobility reveals the 
mechanism underlying slow-rising AMPAR EPSCs and 
the diffusion coefficient in the synaptic cleft. Neuron 
42, 757-771. (doi:10.1016/j.neuron.2004.04.003) 

26. Lehre KP, Rusakov DA. 2002 Asymmetry of glia near 
central synapses favors presynaptically directed 
glutamate escape. Biophys. J. 83, 125-134. 

(doi : 1 0. 1 01 6/S0006-3495 (02) 75 1 54-0) 

27. Takumi Y, Ramirez-Leon V, Laake P, Rinvik E, 
Ottersen OP. 1999 Different modes of expression 
of AMPA and NMDA receptors in hippocampal 
synapses. Nat. Neurosci. 2, 618-624. 
(doi:10.1038/10172) 

28. Nusser Z, Lujan R, Laube G, Roberts JD, Molnar E, 
Somogyi P. 1998 Cell type and pathway 
dependence of synaptic AMPA receptor number 
and variability in the hippocampus. Neuron 21, 
545 - 559. (doi:1 0.1 01 6/S0896-6273 (00)80565-6) 

29. Sabatini BL, Svoboda K. 2000 Analysis of calcium 
channels in single spines using optical flulctuation 
analysis. Nature 408, 589-593. (doi:10.1038/ 
35046076) 

30. Savtchenko LP. 2007 Bilateral processing in chemical 
synapses with electrical 'ephaptic' feedback: a 
theoretical model. Math. Biosci. 207, 113-137. 
(doi:10.1016/j.mbs.2006.09.016) 



31. Wadiche Jl, Jahr CE. 2001 Multivesicular release at 
climbing fiber-Purkinje cell synapses. Neuron 32, 
301-313. (doi:10.1016/S0896-6273(01)00488-3) 

32. Lester RAJ, Tong G, Jahr CE. 1993 Interactions 
between the glycine and glutamate binding sites of 
the NMDA receptor. J. Neurosci. 13, 1088-1096. 

33. Harris KM, Stevens JK. 1989 Dendritic spines of CA 1 
pyramidal cells in the rat hippocampus: serial 
electron microscopy with reference to their 
biophysical characteristics. J. Neurosci. 9, 
2982-2997. 

34. Schikorski T, Stevens CF. 1997 Quantitative 
ultrastructural analysis of hippocampal excitatory 
synapses. J. Neurosci. 17, 5858-5867. 

35. Rusakov DA, Kullmann DM. 1998 Extrasynaptic 
glutamate diffusion in the hippocampus: 
ultrastructural constraints, uptake, and receptor 
activation. J. Neurosci. 18, 3158-3170. 

36. Scimemi A, Fine A, Kullmann DM, Rusakov DA. 2004 
NR2B-containing receptors mediate cross talk 
among hippocampal synapses. J. Neurosci. 24, 
4767-4777. (doi:10.1523/JNEUR0SCI.0364-04.2004) 

37. Oertner TG, Sabatini BL, Nimchinsky EA, Svoboda K. 
2002 Facilitation at single synapses probed with 
optical quantal analysis. Nat. Neurosci. 5, 657-664. 

38. Nimchinsky EA, Yasuda R, Oertner TG, Svoboda K. 
2004 The number of glutamate receptors opened by 
synaptic stimulation in single hippocampal spines. 
J. Neurosci. 24, 2054-2064. (doi:10.1523/ 

J N EU ROSCI .5066-03 .2004) 

39. Christie JM, Jahr CE. 2006 Multivesicular release at 
Schaffer collateral-CA1 hippocampal synapses. 

J. Neurosci. 26, 210-216. (doi:10.1523/JNEUR0SCI. 
4307-05.2006) 

40. Heine M, Groc L, Frischknecht R, Beique JC, Lounis 
B, Rumbaugh G, Huganir RL, Cognet L, Choquet D. 
2008 Surface mobility of postsynaptic AMPARs tunes 
synaptic transmission. Science 320, 201-205. 
(doi:10.1126/science.1152089) 

41. Groc L, Heine M, Cousins SL, Stephenson FA, Lounis 
B, Cognet L, Choquet D. 2006 NMDA receptor 
surface mobility depends on NR2A-2B subunits. 
Proc. Natl Acad. Sci. USA 103, 18 769-18 774. 
(doi:10.1073/pnas.0605238103) 

42. Kennedy MB. 2000 Signal-processing machines at 
the postsynaptic density. Science 290, 750-754. 
(doi:10.1126/science.290.5492.750) 

43. Chen X, Winters C, Azzam R, Li X, Galbraith JA, 
Leapman RD, Reese TS. 2008 Organization of the 
core structure of the postsynaptic density. Proc. Natl 
Acad. Sci. USA 105, 4453-4458. (doi:10.1073/pnas. 
0800897105) 

44. Masugi-Tokita M, Tarusawa E, Watanabe M, Molnar 
E, Fujimoto K, Shigemoto R. 2007 Number and 
density of AMPA receptors in individual synapses in 
the rat cerebellum as revealed by SDS-digested 
freeze-fracture replica labeling. J. Neurosci. 27, 
2135-2144. (doi:10.1523/JNEUR0SCI.2861-06.2007) 

45. Rusakov DA, Kullmann DM, Stewart MG. 1999 
Hippocampal synapses: do they talk to their 
neighbours? Trends Neurosci. 22, 382-388. (doi:10. 
1 01 6/S01 66-2236(99)01 425-3) 



